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Abstract 

In recent years, palladium has proven to be a crucial component for devices ranging from nan- 
otube field effect transistors to advanced hydrogen storage devices. In this work, I examine the 
phonon dispersion of fee Pd using first principle calculations based on density functional perturba- 
tion theory. While several groups in the past have studied the acoustic properties of palladium, this 
is the first study to reproduce the phonon dispersion and associated anomaly with high accuracy 
and no adjustable parameters. 

PACS numbers: 63.20.D, 71. 18,71.20 
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I. INTRODUCTION 



Recent interest in hydrogen storage systems and nanoscale devices has highhghted the 
crucial role palladium plays in a wide range of systems. Hydrogen sensors based on Pd 
nanowiresi show both fast response and low power requirements. Recent experimental work 
indicates that Pd leads provide ohmic contacts for nanotube field effect transistors, a fea- 
ture crucial for large scale device integration^. In addition to nanoscale applications, bulk 
palladium also presents interesting properties that have fascinated researchers for years. It 
possesses a high magnetic susceptibility and sits on the edge of magnetism. Perhaps due to 
this, no measurable superconductivity has been found in the system for any temperature. 

Phonon scattering events play a crucial role in both superconductivity in bulk systems 
and transport in nanoscale interconnects. In this work, I examine the phonon dispersion 
of palladium from a first principles perspective. Recent advances in density functional 
perturbation theory (DFPT) have made it possible to examine the acoustic properties of 
materials at a level of accuracy previously reserved only for electronic properties'*^. Prior to 
the development of this approach, researchers were forced to use a frozen phonon technique 
that required large supercells or phenomenological approaches that relied on numerous fitting 
parameters. 

Interest in the acoustic properties of palladium grew after Miiller and Brockhouse ex- 
amined the phonon dispersion of palladium in detail using inelastic neutron scattering^'^'^. 
They observed a change in the slope of the [CCO] transverse acoustic, TAi, branch of the dis- 
persion curve around g- vector ^[0.35, 0.35, 0] that differed from phonon anomalies observed 
in other metals. The fee lattice constant in this case is given by a. The anomaly extends 
over a broad range of wave vector space and also decreases rapidly with temperature. These 
features initially cast doubt on whether this feature was a Kohn anomaly brought about by 
nesting within parallel sections of the Fermi surface^. 

For small q- vectors in metals, itinerant electrons are able to effectively screen the positive 
ionic charge revealed by lattice vibrations. However once the phonon wavevector, q, spans 
different Fermi sheets and links different electronic states, k2 = ki + q, the electrons are 
no longer able to effectively shield the induced positive charge and the acoustic properties 
undergo a marked change, the Kohn anomaly-. Making use of Fermi surface data from 
augmented plane wave calculations, Miiller found transition vectors between parallel sheets 



2 



of the fifth band hole surface that supported the presence of a Kohn anomaly^ in Pd. Later 
calculations based on the generalized susceptibility, xo il) i^i Pd also indicated that nesting 
between sheets of the fifth band could contribute to the extended anomaly^. This work also 
predicted the presence of a weaker Kohn anomaly at q = ^[0.28, 0.28, 0.28] in the TAi 
branch along the [111] direction. However, Maliszewski et. al. noted that while the sharp 
peak in (q) does occur in the same wave vector range as the [110] phonon anomaly, it 
does not correspond to the behavior of the experimental anomaly in the [110] direction^^. 
They suggested that spin fluctuations or paramagnons could be playing a crucial role in the 
structure of the [111] anomaly. 

Several phenomenological approaches have provided fitting routines to generate the palla- 
dium phonon dispersion. These include studies using a six parameter screened shell model^, 
a three parameter Morse potential approach^, and a four parameter Coulomb potential 
techniquei^ among others. While these approaches can provide general agreement with the 
phonon dispersion of palladium, it should be noted that these techniques do not reproduce 
the phonon anomaly in the [110] direction. The large number of fitting parameters inherit 
in these schemes also limits the predictive ability of these approaches. One study based on 
a pseudopotential approach with a short range pair potential for d-d orbital interactions did 
show evidence for a phonon anomaly in the TAi branch in the [110] directioni^. However, 
the predicted anomaly was shifted to smaller phonon wavevectors and showed a much more 
dramatic phonon softening than the measured anomaly. 

Two recent works have examined acoustic properties of palladium in the framework of 
density functional perturbation theory. Savrasov and Savrasov looked at electron-phonon 
interactions for a range of materials using a full potential linear muffin tin orbital approach 
(FP-LMTO)^S. Their calculations reproduce the general features of the experimental phonon 
dispersion of palladium, however, the phonon wavevector sampling used was too coarse to 
reveal phonon anomalies. A recent estimate of the superconducting transition temperature 
for palladium under pressure also made use of phonon dispersion curves and did not observe 
any phonon anomalies^. 

In this work, I present high resolution calculations for the phonon dispersion of fee pal- 
ladium. I compare direct calculations at q-vectors along high symmetry lines with results 
based on interatomic force constants. In particular, I focus on the phonon anomaly in the 
[110] direction and show that DFPT is able to reproduce this feature with high fidelity. 
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Since this anomaly is due to phonon-electron interactions, I also determine the location of 
the Kohn anomaly based on Fermi surface analysis and the two different techniques show 
good agreement. I have also examined the phonon dispersion along the [111] direction and 
I do not find evidence for the phonon anomaly previously predicted by Freeman et. al.—. 



II. DENSITY FUNCTIONAL PERTURBATION THEORY 

For a periodic lattice of atoms, the only portions of the crystal Hamiltonian that depend 
on the ion positions are the ion-electron interaction terms, Vfj(r), and the ion-ion interaction 
term, ii^jo„(R). The force on a given atom consists of a component that is responding to 
the surrounding electron charge density, n(r), and another term that accounts for ion- ion 
interactions. The interatomic force constants used to determine phonon frequencies in the 
system are found by differentiating the force, F/ on a given ion, /, by the change in position 
of ion J. 



_ -dFj _ fdn{r)dVR{r) 
^^-^^ - W i dRj dRj 

The interatomic force constant matrix (IFC) above consists of seperate electronic and 
ionic contributions. The ionic contribution can be related to an Ewald sum and calculated 
directlyi^. 

There are two factors that are critical for the determination of the electronic contributions 
to the IFC. The first is the ground state electron charge density, ^^^(r). The second is 
the linear response of the ground state electron density to a change in the ion geometry. 
These quantities can be calculated directly within the density functional framework without 
resorting to fitting data from experiment. 

While first principle calculations readily determine the equilibrium electron charge den- 
sity, nji{r), some additional effort is required to determine driRiv) /dRj and the correspond- 
ing IFC. However, the computational effort required is on the same order as the standard 
ground state energy calculation. Interatomic force constants for periodic structures can be 
determined through the use of density functional perturbation theory-i^. In this approach, 
the Kohn-Sham equations for the charge density, self-consistent potential and orbitals are 



linearized with respect to changes in wave function, density, and potential variations. For 
phonon calculations, the perturbation, AVion, is periodic with a wave vector q. This per- 
turbation results in a corresponding change in the electron charge density, An{r). Since the 
perturbation generated by the phonon is periodic with respect to the crystal lattice, we can 
Fourier transform the self-consistent equations for the first order corrections and consider the 
Fourier components, An{c{+ G), of the electron density change, where G is any reciprocal 
lattice vector. In k-space, this problem can now be completely addressed on a sufficiently 
dense k-point grid in the unit cell Brillouin zone and fast Fourier transforms can be used 
to convert back to the real space form. This removes the need to examine frozen phonon 
supercells and greatly accelerates calculations. This approach has been successfully applied 
to reproduce the phonon dispersion of materials ranging from semiconductors, metals, and 
multilayers (see the review by Baroni et. al^ for examples). 

III. RESULTS 

A. Choice of Pd Pseudopotential 

Total energy and phonon calculations were performed on fee palladium using the first 
principles package Quantum-ESPRESSO^*^. Several varieties of local density approximation 
(LDA) and generalized gradient approximation (GGA) pseudopotentials were considered 
in this study including ultrasoft Perdew-Burke-Ernzerhof (PBE) pseudopotentials^^ as well 
as Perdew-Zunger (LDA) ultrasoft pseudopotentials22. In both sets, the effect of including 
semi-core d states was also considered. I found that all GGA pseudopotentials predicted 
a magnetic ground state for fee Pd. In addition, LDA pseudopotentials that did not take 
into account the semi-core d state predicted a magnetic ground state as well. Bulk palla- 
dium is not ferromagnetic, although previous works have indicated a small expansion of the 
crystal can induced a finite magnetic moment^. The LDA pseudopotential that includes 
the semi-core d state provided the physical paramagnetic palladium ground state and is 
used throughout the remainder of the study. The inability for GGA Pd pseudopotentials to 
provide the correct non-magnetic ground state has also been recently found by other groups 
using different first principle approachea^^"^. It is also interesting to note that a recent 
work examining the phonon dispersion in fee palladium found good agreement with LDA 
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pseudopotentials, but poor agreement with GGA pseudopotentials^^. Since the Pd GGA 
pseudopotential overestimates the equihbrium unit cell volume, this leads to an overall shift 
in the predicted phonon dispersion to lower frequencies. 

B. Lattice Dynamics in Palladium 

Total energy calculations made use of a 34 Ry plane wave cut-off and a 348 Ry energy 
cut-off for the electron density. A 16x16x16 shifted Monkhorst-Pack k-point grid^^ was used 
in all calculations and electron occupation was treated with a Methfessel-Paxton smearing 
parameter^^, a, of 0.01 Ry. The smearing parameter plays the role of a fictitious electronic 
temperature and acts to smear out discontinuities found in the density of states at zero 
temperature. The exact result is obtained in the limit as a ^ 0. In most cases, the phonon 
frequencies are fairly insensitive to a, however near Kohn anomalies I found that very small 
a values are often required. 

Phonon frequencies were calculated directly for various phonon wavevectors along high 
symmetry lines using density functional perturbation theory and are shown in Figure [H 
The phonon dispersion was also calculated based on the interatomic force constants obtained 
from a 8x8x8 Monkhorst-Pack g-vector grid. The acoustic sum rule (ASR) was used for the 
interatomic force constant calculations, but this had a fairly small effect of the predicted 
phonon frequencies For example, the maximum difference in phonon frequencies between 
calculations with and without ASR occurred at the P point (Acu = 0.07 THz). In the region 
of the [110] Kohn anomaly, the frequency difference was 0.003 THz or less. The direct 
and ASR interatomic force constant results are shown in comparison to phonon frequencies 
determined in a previous experimental study at 120 K^. The overall agreement for the direct 
phonon calculations is good with an average 5% difference between predicted and measured 
values. The 8x8x8 g-vector IFC results shows good agreement with the direct calculations 
for all of the symmetry lines, but have some difficulty reproducing experimental results 
at low frequencies in the acoustic branches. I found coarser g-vector grids were unable to 
reproduce the experimental phonon dispersion curves. 

It should be stressed that the phonon dispersion curves from both the direct calculations 
and interatomic force constants are calculated without any adjustable parameters and make 
no reference to experimental data. There are two clear sources for differences between the 
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measured and calculated phonon frequencies. The first is due to the fact that experimental 
studies were performed at 120 K, while first principle calculations assume zero temperature. 
Thermal expansion of the unit cell will act to shift the phonon frequencies lower. In ad- 
dition, the local density functional approximation (LDA) has a tendency to overbind and 
predict slightly smaller lattice constants than those observed experimentally. Using LDA, 
the equilibrium lattice constant for fee Pd was determined to be 3.88 A, while the measured 
lattice constant is slightly larger, 3.89 A. The smaller predicted lattice constant will lead to 
a more compressed unit cell and a general shift to higher phonon frequencies. This trend of 
overestimating the phonon frequency is greatest in the longitudinal acoustic branch. 

While there are some slight differences between the current work and previous experi- 
mental studies, it is important to note that the direct and IFC calculations replicate the 
phonon anomaly observed in the [110] direction. Previous parameter based approaches have 
been unable to capture this anomaly. This feature in the TAi branch is shown in Figure [2] in 
comparison with experimental data at 8 The slope of the calculated phonon dispersion 
changes at approximately q = ^[0.30, 0.30, 0], indicating a Kohn anomaly. The slope of 
the experimental phonon dispersion changes at a slightly higher q- vector, q = ^[0.35, 0.35, 
0] . This small difference in the calculated and experimental q- vector for the Kohn anomaly 
could be related to the smaller lattice constant predicted by density functional techniques. 

As an additional check to determine whether this feature is a Kohn anomaly, I examined 
the TAi branch as a function of the smearing parameter, a. Since this parameter serves as 
a fictitious temperature, increasing it will affect virtual transitions between occupied and 
empty states in an energy range a around the Fermi energy^. This will govern how effectively 
electrons are able to shield motion in the ion lattice and should lead to the disappearance 
of the anomaly. The TAi branch for smearing parameters ranging from 0.01 Ry to 0.10 Ry 
are shown in Figure [3] in comparison with the experimental values. It is clear that as the 
smearing parameter increases, the Kohn anomaly disappears. For phonon wavevectors far 
from the Kohn anomaly, like q = ^[0.70, 0.70, 0], adjusting the smearing parameter has 
little effect on the TAi phonon frequency. 

The Fermi surface for fee Pd was extracted from the calculation as well. There are three 
bands that contribute at the Fermi energy (Fig. 111(a)). The fourth band barely overlaps 
with the Fermi energy and only exhibits small cusps at the center of the [100] faces. The 
sixth band is centered in the Brillouin zone and although faceted, has high s character and 



7 



is similar to a free electron Fermi surface. The Fermi surface for the li-like 5th band in fee 
palladium is shown in Figure 111(b). 

The fifth band Fermi surface is dominated by wide sheets that sit in the [100] faces of the 
Brillouin zones. These parallel sheets provide important contributions to the Kohn anomaly 
in the [110] direction. This can be seen by examining the intersection of the fifth band with 
the kz = plane in the central Brillouin zone and neighboring Brillouin zones (Fig. [5]). 
The fifth band Fermi surface forms squares in the kz = plane that can be linked with 
q = ^[0.30, 0.30, 0] nesting vectors. These nesting vectors occur at the phonon wavevector 
calculated for the Kohn anomaly in the [110] direction. 

Based on analysis of the Fermi surface, I found that the s character sixth band also 
provides contributions to the Kohn anomaly in the [110] direction. As shown in Figure El the 
sixth band Fermi surface has fairly wide flat valleys in the [110] directions that are essential 
for nesting vectors and Kohn anomalies. The portions of the sixth band that intersect the 
kz = plane are marked with blue lines for the center Brillouin zone and neighboring ones 
in Figure [71 In the [110] directions, these cross-sections have almost no curvature, making 
them ideal for Fermi nesting. A representative nesting vector based on the Fermi surface 
between adjacent Brillouin zones is denoted by a black arrow. The average nesting vector 
based on the Fermi surface in the [kx,ky] plane is given by q = ^[0.30, 0.30, 0] and is 
in good agreement with the Kohn anomaly calculated from density functional perturbation 
theory. The observed Kohn anomaly is fairly broad due to the fact that parallel segments 
of the Fermi sheets off the [kx, ky] plane also contribute to the anomaly at slightly different 
q vectors. 

The phonon frequencies along the [111] direction were calculated directly for both the 
transverse and longitudinal branches in an effort to resolve the phonon anomaly predicted 
by Freeman et al^. However, in this case, there was no clear evidence of a phonon anomaly 
in either acoustic branch. 

IV. CONCLUSION 

In this work, I have considered the lattice dynamics of fee palladium based on density 
functional perturbation theory. Agreement between direct phonon calculations along sym- 
metry lines with current experimental data for the full phonon spectrum of palladium is 



8 



good. In particular, I was able to reproduce the extended phonon anomaly in the [110] di- 
rection that had been observed experimentally. In addition, I found no evidence for phonon 
anomalies in other branches of the phonon dispersion. I also examined the phonon dispersion 
derived from interatomic force constants and looked at the influence of g-vector sampling. 
The interatomic force constant approach with a 8x8x8 g-vector grid was able to provide 
good agreement for general features of the palladium phonon dispersion. 
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List of Captions 

Figure 1: (color online) Phonon dispersion for fee palladium along different symmetry 
lines. The phonon dispersion based on the 8x8x8 q-vector grid interatomic force con- 
stants is given by black lines. The results from the direct phonon calculations are 
denoted by black triangles and shown in comparison with experimental results from 
Miiller and Brockhouse^ taken at 120 K (red circles). 

Figure 2: (color online) Phonon dispersion for the TAi branch of fee Pd along the [CCO] 
line from F to X. The direct calculations (red triangles) and IFC-based results (blue 
line) are shown in comparison with the experimental results from Miiller- taken at 8 
K (black circles). The location of the phonon anomaly is denoted by the arrow on the 
graph. 

Figure 3: (color online) Phonon dispersion for the TAi branch of fee Pd is shown along 
the [CCO] line from F to X for different values of the smearing parameter. The results 
from calculations using cr = 0.01 Ry (open diamonds), cr = 0.05 Ry (black squares), 
and (J = 0.10 Ry (red circles) are shown in comparison with the experimental results 
from Miiller- taken at 8 K (triangles). 

Figure 4: (color online) The Fermi surface as viewed head on from the [100] direction is 
shown in (a). The Fermi surface for the fourth band forms small cusps in the [100] 
faces and is shown in dark blue. The Fermi surface for the fifth band is denoted by 
a purple interior and green exterior and forms petals center on the [100] faces. The 
sixth band with high s character is located at the center of the Brillouin zone and is 
shown in light blue. In part (b), the Fermi surface for the fifth band is shown alone 
for greater clarity. 

Figure 5: (color online) The intersection of the fifth band Fermi surface (orange line) 
through the kz = plane is shown. The Fermi surface intersections in relevant neigh- 
boring Brillouin zones are also shown (blue lines). The nesting vectors, q = ^[0.30, 
0.30, 0] between adjacent sections of the Fermi surface is denoted by black arrows. 

Figure 6: (color online) The Fermi surface for the sixth band is shown in light blue within 
the fee Brillouin zone (red lines). This band has a high s character and peaks in the 
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[Ill] and [100] directions. A dashed line indicates the center of the wide flat vaUeys in 
the [110] directions. These regions are crucial for the nesting vectors which generate 
the Kohn anomaly in the [110] direction. 

Figure 7: (color onhne) A cross-section of the sixth band Fermi surface (blue lines) through 
the kz — plane is shown. The Fermi surface cross-sections from neighboring Brillouin 
zones are also shown. The nesting vector, q = ^[0.30, 0.30, 0] between adjacent 
sections of the Fermi surface is denoted by a black arrow. 
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